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I. INTRODUCTION 

In a recent paper [lj , the 1 + 3 formalism [I! has 
been used to develop a general framework for studying 
axially symmetric dissipative fluids. More recently, the 
above mentioned framework has been applied to study 
the perfect and geodesic, fluid case without imposing 
from the beginning the shear-free condition. However, as 
a result of this study, it appears that all possible models, 
compatible with the assumed line element, are in fact 
shear-free, and according to the result obtained in 0, 
they are Friedman-Robertson-Walker (FRW), i.e. they 
do not produce gravitational radiation. 

We recall that we define a state of intrinsic gravita¬ 
tional radiation (at any given point), to be one in which 
the super-Poynting vector does not vanish for any unit 
timelike vector |8j-jl0]. Then since the vanishing of the 
magnetic part of the Weyl tensor implies the vanishing of 
the super-Poynting vector, it is clear that FRW does not 
produce gravitational radiation. It is also worth recall¬ 
ing that the tight link between the super-Poynting vector 
and the existence of a state of radiation, is firmly sup¬ 
ported by the relationship between the former and the 
Bondi news function [H], [lj[ (see [h| for a discussion on 
this point). 

The fact that a non-dissipative source, does not radi¬ 
ate gravitational waves, becomes intelligible if we recall 
that radiation is an irreversible process. This emerges 
at once if absorption is taken into account and/or Som- 
merfeld type conditions, which eliminate inward traveling 
waves, are imposed 0 - Accordingly, an entropy gener¬ 
ator factor should be present in the description of the 
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source, a factor which of course is absent in a perfect 
fluid, and even more so in a collisionless dust. In other 
words, we should expect the irreversibility of the process 
of emission of gravitational waves, to be represented in 
the equation of state through an entropy increasing (dis¬ 
sipative) factor. 

From the comments above, it becomes clear that the 
simplest fluid distribution which we might believe to be 
compatible with a gravitational radiation, is a dissipative 
dust under the geodesic condition. Such a case will be 
analyzed in this work. 

We shall consider separately the two possible subcases, 
namely: the fluid distribution is assumed, from the be¬ 
ginning, to be vorticity-free, or not. 

In the former case, it is shown that the vanishing 
vorticity implies the vanishing of the heat flux vector, 
and therefore, as shown in (6], the resulting spacetime is 
FRW. 

In the latter case, it is shown that the enforcement 
of the regularity conditions at the center, implies the 
vanishing of the dissipative flux, leading also to a FRW 
spacetime. 

Thus all possible models, sourced by a dissipative 
geodesic dust fluid, belonging to the family of the line ele¬ 
ment considered here, do not radiate gravitational waves 
during their evolution, unless regularity conditions at the 
center of the distribution are relaxed. 

We shall discuss about this result and its relationship 
with the peculiar evolution of the vorticity under the 
geodesic condition. 

This work heavily relies on the general framework de¬ 
veloped in [l|, keeping the same notation, and just re¬ 
ducing the general equations to the particular case con¬ 
sidered here. These will be presented in an Appendix. 
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II. THE DISSIPATIVE, GEODESIC, DUST 

Let us consider axially and reflection symmetric 
geodesic dust distributions (not necessarily bounded). 
For such a system the line element in “Weyl spherical 
coordinates”, is assumed as in Q: 

ds 2 = -dt 2 + B 2 (dr 2 + r 2 d0 2 ) + 2 G(r, 6)dtd8 + C 2 d<t> 2 , 

(!) 

where B and C, are positive functions of t, r and 9 , 
and G depends only on r, 6. We number the coordinates 
x° = t 7 x 1 = r, x 2 = 6 , a; 3 = <j>. 

The energy momentum tensor in the “canonical” form 
reads: 

T a f) = lBV a Vp + q a Vp + qpV a , ( 2 ) 

where as usual, //, q a and Vp denote the energy density, 
the heat flux vector and the four velocity, respectively. 

The heat flux vector can be decomposed, in terms of 
two scalar functions qi,qu, as (see [l| for details) 

q^i = qiK^ ( 3 ) 

with vectors K and L having components: 

K a = (0, B, 0,0); L a = (0, 0, V B 2 r 2 + G 2 , 0). (4) 
Next, the shear tensor 

O’ci/3 = +(a;/3) + a (a.Vp) — (5) 

where 

h a p =6$ + V a Vp, ( 6 ) 


whereas the vorticity may be described, either by the 
vorticity vector w“, or the vorticity tensor defined 
as: 

= \ rtaPv, V^ V v = ^ n a p» v v v , ( 11 ) 

where Q a p = V\ a -p-\ +a\ a Vp\, and r} a p^ v denotes the Levi- 
Civita tensor. 

We find a single component different from zero, pro¬ 
ducing: 

fl a p = Q(L a Kp — LpK a ), (12) 

and 

U> a = —ClS a . (13) 

with the scalar function Q given by 

n =- , f - ■ (14) 

2B\J B 2 r 2 + G 2 

Observe that from m and regularity conditions at 
the centre, it follows that: G = 0«-!1=0. 

Regularity conditions at the center, necessary to ensure 
elementary flatness in the vicinity of the axis of symmetry 
require that as r s=s 0 

n = ^Q<")(t,0)r n , (15) 

n> 1 

implying, because of (fill) , that in the neighborhood of 
the center 


may be defined through two scalar functions, as: 


& a p — ^ (2(7/ + (77/) (K a Kp ^ h a( g) 
+ “(2(777 + &i)(L a Lp — —h a p). 


(7) 


The above scalars may be written in terms of the metric 
functions and their derivatives as (see 0 ): 


a i = 


B 2 r 2 + 2 G 2 B C 
B 2 r 2 + G 2 B C 


( 8 ) 


(77 - (77/ = 


3G 2 B 
B 2 r 2 +G 2 B’ 


G=^& n \6)r n , 


(16) 


n> 3 


This last result in turn implies that as r aproaches 0, 


(77 - (7// = [°/ n) (M) - r n . (17) 


i>4 


Now, from the elementary flatness condition, we may 
write, as r —> 0, 


c ss n{t,9), 


(18) 


(9) implying 


For the other kinematical variables (the expansion and 
the vorticity) we have: 


2 B 2 r 2 + G 2 B ( C 
B 2 r 2 +G 2 B + C’ 


( 10 ) 


C' ~ j(t, 9), C g ss (19) 

where q(f, (9) is an arbitrary function of its arguments. 

We shall treat separately, the two possible cases, with 
and without vorticity. 

Let us first consider the restricted case G = 0 4=> f2 = 0. 
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III. THE VORTICITY FREE CASE 


In order to prove the main result of this work, i.e., 
that a source represented by a geodesic, dissipative, dust, 
endowed with axially and reflection symmetry, is just 
Friedman-Robertson-Walker (FRW), we shall first con¬ 
sider the case when the dust is irrotational. 

Thus, if we assume the vorticity to vanish, then, from, 
m El ESI) we get 


0 = 2 § + §> CT/ = CT// = ^ = ;§ - §’ ( 20 ) 


implying: 


. „„ 3C 3 B _ 

e = 2<, + — = — -„. 


( 21 ) 


In this case (1A.1I) produces the following two equations 


1 

Jr? 



,e — o,e — cr- 


3C ( 


C 


8irqn, 


(31) 


<r /era C,e\ _ (5-C),0 

2 Br\ 5 + C ) ~ 2 BrC ’ 


(32) 


H 2 = 


(7 

2B 



(dcy 
2 BC ' 


Combining (l32l) with (1331) . we obtain 

(HiBrCy + {H 2 BC),e = 0, 


(33) 


(34) 


where Hi , H 2 are the two scalars defining the magnetic 
part of the Weyl tensor (see Eq.(37) in [1]). 

From (13U1) and (13TT) . it follows that 


Qi + ^-(40 + cr) — 0, (22) 


2(0 + a)' = QeSlL + 24tt qi B, (35) 

o 


qn H——(40 + ct) — 0. (23) 

Using in (1771 [Til) , we obtain at once 

qiB 3 rC = ai(r,6)] qnB 3 C = a 2 (r,0), (24) 

implying 

qi = q{r,9)q n , (25) 

where oi(r, 9 ), a 2 (r, 9) are arbitrary functions of their ar¬ 
guments, and 


2(0 + a) t g = _|_ 24irq n Br, 

which by virtue of m, become 


B\ B C’ C' „ „ 

b) = bc'c +m = 0 ’ 


B\ B Co Co 
B =B^-^ + 87rqiirB ’ 


(36) 


(37) 


(38) 


q(r,9) 


Oil 

a 2 r 


(26) 


Next from (IA.6I) in the vorticity free case, 
(qjBrC)' + (quBC), g = 0. 
The combination of (1771) with (1271) produces 
B' B t e a'l + a 2 ,e 

— Oil H-— 02 — ---, 


we obtain: 

(27) 

(28) 


where (1771) has been used. 

Feeding back these two last equations into m EH)) 
and using (1771) . we obtain 


Hi = A-nqn, H 2 = -4 7r<?/. (39) 

Next, combining (1A.111lAjTl) with (13711551) we obtain 
(qnBrCy - { qi BC), e = 0. (40) 

From 433 ), 433 , using © and (1771) . we obtain 


which implies 


B = T(t)B(r,9), (29) 


TC - TC' 
TC,e - TC, e 


q{r, 9) 


q{r,9 ) 


(41) 


where T and B are arbitrary functions of their argu¬ 
ments. 

Next, (IA.71 IA.81 IA.9I [OH . read, in this case, respec¬ 
tively: 

^20' -a' - = 8nqj, (30) 


whose first integral may be written as 

T(t)R(r,6) = C' -C, e q(r,9), (42) 

where R{r, 9) is an arbitrary integration function. 

The above equation, in turn, produces 

C = T(t)C(r,6), (43) 
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with 


R(r,9)=C'-C, e q(r,9). (44) 

Feeding back (S3) into ill) and (l38l) . it follows at once 
that qj = qii = 0, which by virtue of (1551) implies that the 
magnetic part of the Weyl tensor vanishes. Furthermore, 
since the fluid is geodesic and non-dissipative it must be 
FRW, as shown in 

We shall next consider the case where the vorticity is 
not assumed to vanish from the beginning. 


IV. THE NON-VANISHING VORTICITY CASE 


In the previous section we have shown that, if the van¬ 
ishing vorticity condition is assumed ab initio , then, the 
fluid configuration is just FRW. We shall now prove that 
this is the case even if such a condition is not imposed 
from the beginning. 

Thus, let us now relax the vorticity-free condition. 

In this case we shall proceed as follows. Contracting 
(1A. 141) with vectors K and L, we obtain respectively 


£lqn + 


nT r 

Bt 

qi°i 


1 1 
- + -D t 
T 2 




-6 e > 9 ' 


= 0 , 


(45) 


and 


- SI qj + 


nL a T c 

T 

quern 


i + I A 

T 2 




= 0 , 


-6 e i 9 " 

(46) 


with D t f = f, p VP. 

Let us now impose the regularity conditions at the cen¬ 
ter in (l45l l46|l . 

Close to the center, we may write 


qiI n= ]T q%\t,9)rt n Ht,9)r n , (47) 

n>N 


where N > 1, and which implies because of m 

m= E 0)r n . (48) 

n>N—l 


Feeding back m into P5]) . it follows at once that as 
r « 0, 


^=E^ n) M) r "’ ( 49 ) 

n>N 

where the fact has beeen used that T r and qi have the 
same behaviour at the center, and of course, all physical 
and geometrical variables are regular at r ss 0. 


But then, using (l49l) in (1461) it follows that 

qn = E 9/’r ) (*i 6 ') r "> (50) 

n>./V+l 

which of course is more restrictive than (l48|) . 

We must then assume m instead of (1551) . in which 
case (1451) produces 

qi = E Qi n) (t,6)r n , (51) 

n>N+2 

implying because of (f46|) 

qn = E 1u\t,d)r n . (52) 

n>JV+3 

Going through this cycle, once and again, as many 
times as desired, it is obvious that 

Qi = E qn = E 1ii\t,S)r n , (53) 

n>N n>N 


where now N is an arbitrarily large number. In other 
words, not only <?/.//, but also all their r-derivatives of 
any order, vanish at the center. This in turn implies that 
we can proceed exactly as in subsections II.B and II.C 
in Q, and prove that close to the center we also have 
H[ n) = H^ n) = 4 n) = 4? = = °n = tt (n) = 0 for 

any value of n > 0. 

Then, assuming that qiji are of class C u , we can an¬ 
alytically continue the zero value at the center to the 
whole configuration, i.e. we have a non-dissipative fluid. 
In this case as we know from Q, the system is a FRW. 


V. CONCLUSIONS 

At the light of the results obtained in Q, we started 
this research harboring the hope, that the inclusion of the 
dissipative term in the geodesic fluid distribution, would 
be enough to provide us with a well behaved source of 
gravitational radiation. However, the final result indi¬ 
cates that the geodesic condition is too restrictive, lead¬ 
ing to the impossibility of ensuring regularity conditions 
at the center, in the presence of a dissipative flux (i.e. 
bringing us back to the case analyzed in fH). 

But, could we guess such a result, before carrying out 
the discussion presented here, only based on the geodesic 
condition? As we shall see now, the answer to this ques¬ 
tion, might be yes. 

Indeed, in the fl = 0 subcase, the above mentioned 
result might have been inferred based on the link between 
vorticity and gravitational radiation E3- 

To see how this affirmative answer might have been 
guessed, in the second subcase (i.e. when the vorticity is 
not assumed to vanish from the beginning), let us recall 
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the evolution equation for the vorticity scalar (Eq. B5 in 
JT| . for the geodesic case), which reads 

0 + \ (20 + crj + a n )n = 0. (54) 

o 

Now, let us ignore, for a moment, the result obtained 
in this paper, about the fact that the geodesic condi¬ 
tion leads to a FRW spacetime, implying the vanishing 
of vorticity, and let us focus on the equation above. 

From (1541) we see that the vanishing of vorticity at any 
given time, implies its vanishing at any other time after¬ 
wards, or, alternatively, the presence of a non—vanishing 
vorticity in a geodesic fluid distribution, at any given 
time, implies that such a vorticity was there always , 
before that time. In other words, any possible non¬ 
vanishing vorticity in such a system, would be a primor¬ 
dial one. Such a conclusion, of course, is independent on 
whether or not the system is dissipative. 

Now, the situation depicted above prevents the possi¬ 
bility to have an initially static system, which at some 
given moment departs from equilibrium and starts to 
evolve with a non vanishing vorticity, during a finite pe¬ 
riod of time. But then, and based on the link between 
vorticity and gravitational radiation, mentioned before, 
we might conclude, just from (1541) . that gravitational ra¬ 
diation from a physically meaningful system, i.e. one 
which radiates for a finite period of time, is not to be 
expected under the geodesic condition, even though the 
dissipation is not excluded ab initio. 

Of course, all these “hand waving” arguments, as 
sound as they may be, only provide the basis for a con¬ 
jecture, and cannot replace the formal proof presented in 
sections III and IV. 

The following comments are in order before ending: 

1. In J6j, it was shown that, only using the evolu¬ 
tion equation for the vorticity, for the perfect (i.e. 
in thermodynamic equilibrium) and, non -geodesic, 
fluid, the vorticity is also a primordial one. Indeed, 
in this latter case, the evolution equation for the 
vorticity becomes 

n tS V 5 + ^(20 + ai + a H + WT iM )D = 0. (55) 

o 

where T = InT, and T denotes the temperature. 
Thus, even if the fluid is not geodesic, but is non- 
dissipative, the situation is the same as in the 
geodesic case, i.e. the vanishing of vorticity at any 
given time implies its vanishing for any time in the 
future. 

We then noticed that such a result was in agree¬ 
ment with earlier works indicating that vorticity is 
produced by dissipation HU - !!!]’ which was absent 
in that discusssion. Here we do have a dissipative 
term, but only before enforcing the regularity con¬ 
ditions, and nevertheless the nature of the vorticity 
is of the same kind as that of the dissipationless 


case. Of course, once regularity conditions are en¬ 
forced, the resulting fluid is just FRW, implying 
vorticity and dissipation, vanish. 

2. We must recall that the line element (ID, for which 
our results were obtained, is not the most general, 
compatible with the axial symmetry. In particu¬ 
lar, in the spherically symmetric limit, it describes 
only shear-free fluids. This is so because the shear, 
is now sourced by the magnetic parts of the Weyl 
tensor. 

Indeed, from (1551) and (1551) we find that if the mag¬ 
netic part of the Weyl tensor vanishes, then 


aC = ijj(t ), (56) 

where ipit) is an arbitrary integration function. 
The above equation implies, because of the regu¬ 
larity condition (fl8l) . that a = 0. 

3. Models compatible with the geodesic condition do 
exist, if regularity conditions at the center are re¬ 
laxed. 

4. If regularity conditions at the center are main¬ 
tained, then, physically acceptable models require 
the inclusion of, both, dissipative and anisotropic 
stresses terms, i.e. the geodesic condition must be 
abandoned. In this case, purely analytical meth¬ 
ods are unlikely to be sufficient to arrive at a full 
description of the source, and one has to resort to 
numerical methods. 
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Appendix: Summary of equations for the geodesic 
dissipative case 

Below, we shall write only the equations required for 
our discussion, and specialized for the geodesic and dis¬ 
sipative dust, from the framework developed in [l]. 

The equation (A7) in [l[, which comes from the con¬ 
servation law To. a = 0, leads in our case to: 

+ 0a/3 + = 0- (A.l) 

which is the “generalized Euler” equation. 
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From B.10, B.12, B.13, in pj, which are some of the Bianchi identities with the tetrad vectors, we obtain re¬ 
scalar equations resulting from the projections of the spectively 


\&i + ~ &£kl + - 77 - (30 + 077 c 1 ) ~E~ (2 an + 07 ) 

00 y y 


1 


\jB 2 r 2 + G 2 


Hie + H, 


C 


Ei 

B 


C_ _ 2 {Br){Br)' + GG' 
~C ~ 


(' B 2 r 2 + G 2 


47T 

~B 


<h 


4 nqn 


\J B 2 r 2 + G 2 


gE + Ea 

B B 


(A.2) 


7^77 + + VLSkl + (30 + oy — an) + (2oy + an) 


+1 ( H ' 2 + hS\+ Hl _ 

B V C ) yjB 2 r 2 + G 2 


G B \B C 


= —2tt 


q 1 ( B 2 r 2 -(- G 2 )' 4ir(Gqn + gu,e) 


B B 2 r 2 + G 2 


V B 2 r 2 


G 2 


(A.3) 


— 2 (£7 + £ 77 )' - ^(£7 + £ 77 )© — -g-/7( a 1 + <7//) — (2a n + aj) -(2 07 + an) 


1 


VbV + g 2 


H\ g + \ — — < H' 2 + H 2 


(. B 2 r 2 + G 2 )' 


2(B 2 r 2 + G 2 ) 


47TQ/ C' 

'—~c 


Anqn 


2 r 2 
C C, 


G 


\JB 2 r 2 + G 2 \ C C 


A combination of (IA.3I) and (IA.4I) produces 


(A.4) 


~ \&i - nrl^r + ?L£kl - 77 -(30 - oy + a n ) - ^-(oy + 2an) 

00 y y 


+ 


VB 2 


G 2 


H he + H 1 


Ei-r E-A 

c g \b c 


El 

B 


C' (S 2 r 2 + G 2 )' 


G 2(B 2 r 2 + G 2 ) 


47T 

= “S ? 7 


G' (B 2 r 2 + G 2 )' 


G 2{B 2 r 2 + G 2 ) 
whereas from (IA.2I) with (IA.5I) we obtain 


47T 


V'BV 2 + G 2 


Gg// + qn,e + qn f G— + - ■pj- 


(A.5) 


G 


\J~B^r 


2„2 


G 2 


H 2 G' 


+ Hi 


B C 


2B\JB 2 r 2 + G 2 ' \ B C 


4tt I , 

= ^ 97 + 97 


G' (B 2 r 2 + G 2 )' 


G 2(B 2 r 2 + G 2 ) 


47T 


a/b 2 


G 2 


Qn 


G 


C B 
C + B 


Ea 

c 


Ea 

B 


+ Gqn + qn,e 


(A.6) 


Next, from B.6, B.7, B.8, and B.9 of [I], which are obtained by contracting the Ricci identities with the tetrad 
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vectors, we obtain 


\jB 2 r 2 + G 2 




(GG C,b Y 

, 1 I 

{ c + c ) 

+ 3 B | 

\ / . 



+ 0 / 


B 

Be Cg ~ 

~b~~c +g 


i / a-" 

17' + ] + 


-0/7 


C' 

(7 


26" (SV+G 2 ) 

<? + 2 


' B 2 r 2 + G 


( B 2 r 2 + G 2 



WB 2 r 2 + G 2 


/b c V 


(5-cjJ 

- C/7 


|(20 — an) g + G (20 — an)' 
= 8nqn, 


Be 2C e ~ B 2G 

—+ —+ G - + — 


(A.7) 


(A.8) 


^ = -21 ff - n 


G' 


GG' 


G 2(B 2 r 2 + G 2 ) 


+0/ 


H£j.9£-rr?.-<L\ 

B + C G ^B c' 


- on 


6VB 2 r 2 + G 2 


B, e 2C fi ~2B 2 C, 
B G + B G ; 


(A.9) 


H 2 = <1 (oj + 2(7//)' + (7/ 


2G' (Br)(Br)' 

~C ~ B 2 r 2 + G 2 

1 


+ on 


(y_ 2 (Br)(Br)' + 3 GG' 

~C + 


2\J B 2 r 2 + G 2 


£1 /Q — £7 


(^B 2 r 2 + G 2 


4 HI 4 


and from B.ll and B.16 of [l|, which are scalar equations, also obtained from the Bianchi identities 


2 £kl + £kl (20 — oi — an) + — (£/ — £n) + — { H[ + H± 


H 2 ,e + Hi 


\jB 2 r 2 + G 2 


2 C' 2(Br)(Br)' + GG 1 
G 2 (B 2 r 2 + G 2 ) 

2G, e B,* ~(B G N 

G B G [b C 


2t r 


v^bvTg 2 


9/ I G— + — I - Gqj - q/,g 


27T 

IT 


-9// + Qii 


(B 2 r 2 + G 2 )' 


1 


= [87t^ - (£7 + £//)] fl - 


Vbv + g 2 

47T 


"farL {on - a T ) - — \h[ + Hx 


H 2 ,e + GH 2 + H 2 


2{B 2 r 2 + G 2 ) 

2G' , (B 2 r 2 + G 2 )' 

~C + 2{B 2 r 2 + G 2 ) 


a/BV 2 + G 2 


, b , 9 \ 4 tt 

q/,e + 97 - 5 - + 


B 


B 



(A.10) 


(A.ll) 


(A.12) 


































































(A.13) 


Finally, the transport equation for the geodesic case, (Eq.(57) in 0 ) reads: 

i / T ycx\ 

ThXpVt + q^ = ~kW v T v - -kT 2 (^J q», 

then, from a combination of (IA.1 1 ) and (IA. 131) (Eq.(60) in |lj), we obtain in our case. 


— (Oa /3 + ^ap)q^ + —Vo ,T + | — + — D t ln(- 7 ^ 2 ) — g0 \ Qa —0. 


(A. 14) 


[1] L. Herrera, A. Di Prisco, J. Ibanez and J. Ospino, Phys. 
Rev. D 89 , 084034, (2014). 

[2] G. F. R. Ellis Relativistic Cosmology in: Proceedings 
of the International School of Physics “ Enrico Fermi”, 
Course 47: General Relativity and Cosmology. Ed. R. K. 
Sachs (Academic Press, New York and London) (1971). 

[3] G. F. R. Ellis and H. van Ellst, gr-qc/9812046v4 (1998). 

[4] G. F. R. Ellis Gen. Rel. Grav. 41, 581 (2009). 

[5] G. F. R. Ellis, R. Maartens and M. A. H. MacCallum, 
Relativistic Cosmology (Cambridge U. P., Cambridge) 
( 2012 ). 

[6] L. Herrera, A. Di Prisco, J. Ospino and J. Carot Phys. 
Rev. D 91, 024010, (2015). 

[7] L. Herrera, A. Di Prisco and J. Ospino, Phys. Rev. D 
89,127502, (2014). 

[8] L. Bel, Cah. de Phys. 16 , 59 (1962). 

[9] A. Garcla-Parrado Gomez Lobo, Class. Quantum Grav. 
25 , 015006 (2008). 

[10] A Di Prisco, L. Herrera, M. A. H. MacCallum and N.O. 
Santos, Phys. Rev. D 80 , 064031 (2009). 

[11] H. Bondi, M. G. J. van der Burg and A. W. K. Metzner, 
Proc. Roy.Soc. A 269 , 21 (1962). 


[12] R. Sachs, Proc. Roy.Soc. A 270 , 103 (1962). 

[13] L. Herrera, W. Barreto, J. Carot and A. Di Prisco, Class. 
Quantum Grav., 24 , 2645, (2007). 

[14] H. Stephani, D. Kramer, M. MacCallum, C. Honselaers, 
and E. Herlt, Exact Solutions to Einsteins Field Equa¬ 
tions (Cambridge University Press, Cambridge, Eng¬ 
land), (2003), 2nd Ed. 

[15] J. Carot, Class. Quantum Grav. 17 , 2675 (2000). 

[16] G. T. Carlson, Jr. and J. L. Safko, Ann. Phys. (N. Y.) 
128 , 131 (1980). 

[17] L. Herrera, Gen. Rel. Grav. 46, 1654 (2014). 

[18] L. Crocco and Z. Agnew, Math. Mech. 17, 1 (1937). 

[19] E. R. Harrison, Mon. Not. R. Astr. Soc. 147 , 279 (1970). 

[20] A. J. Christopherson, K. A. Malik and D. R. Matravers, 
Phys. Rev. D 79 , 123523 (2009). 

[21] F. Del Sordo and A. Brandenburg, Astr. Astrophys. 528 , 
A145 (2011). 

[22] A. J. Christopherson and K. A. Malik, Class. Quantum 
Grav. 28 , 114004 (2011). 

[23] F. Dosopoulou, F. Del Sordo, C. G. Tsagas and A. Bran¬ 
denburg, Phys. Rev. D 85 , 063514 (2012). 











